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ABSTRACT 

We present the second in a series of papers discussing the thermodynamic properties 
of M 87 and the central regions of the Virgo Cluster in unprecedented detail. Using 
a deep Chandra exposure (574 ks), we present high-resolution thermodynamic maps 
created from the spectra of ^^16,000 independent regions, each with ~1,000 net counts. 
The excellent spatial resolution of the thermodynamic maps reveals the dramatic 
and complex temperature, pressure, entropy and metallicity structure of the system. 
The 'X-ray arms', driven outward from M87 by the central AGN, are prominent 
in the brightness, temperature, and entropy maps. Excluding the 'X-ray arms', the 
diffuse cluster gas at a given radius is strikingly isothermal. This suggests either that 
the ambient cluster gas, beyond the arms, remains relatively undisturbed by AGN 
uplift, or that conduction in the intracluster medium (ICM) is efhcient along azimuthal 
directions, as expected under action of the heat-flux driven buoyancy instability (HBI). 
We confirm the presence of a thick 40 arcsec or ^ 3 kpc) ring of high pressure 
gas at a radius of 180 arcsec (~ 14 kpc) from the central AGN. We verify that 
this feature is associated with a classical shock front, with an average Mach number 
M = 1.25. Another, younger shock-like feature is observed at a radius of ~ 40 arcsec 
(~ 3 kpc) surrounding the central AGN, with an estimated Mach number M > 1.2. 
As shown previously, if repeated shocks occur every ~ 10 Myrs, as suggested by these 
observations, then AGN driven weak shocks could produce enough energy to offset 
radiative cooling of the ICM. A high significance enhancement of Fe abundance is 
observed at radii 350 — 400 arcsec (27 — 31 kpc). This ridge is likely formed in the 
wake of the rising bubbles filled with radio-emitting plasma that drag cool, metal-rich 
gas out of the central galaxy. We estimate that at least ~ 1.0 x 10^ solar masses of Fe 
has been lifted and deposited at a radius of 350 — 400 arcsec; approximately the same 
mass of Fe is measured in the X-ray bright arms, suggesting that a single generation 
of buoyant radio bubbles may be responsible for the observed Fe excess at 350 — 400 
arcsec. 

Key words: X-rays: galaxies: clusters - galaxies: individual: M87 - galaxies: inter- 
galactic medium - cooling fiows 



1 INTRODUCTION 

The Virgo Cluster is the second brightest, extragalactic, ex- 
tended X-ray source in the 0.1 — 2.4 keV ROSAT band, and 
the closest galaxy cluster (16.1 Mpc; Tonry et al. 2001). The 
central galaxy of the Virgo Cluster, M 87, hosts an active 
galactic nucleus (AGN) that exhibits compelling evidence 
for complex interactions with the surrounding cluster gas 



{e.g. Bohringer et al. 1995; Young et al. 2002; Forman et 
al. 2005, 2007). M87 is, therefore, an excellent object with 
which to study AGN driven processes {i.e. AGN feedback) 
that affect the centers of galaxies and galaxy clusters. 

The relatively cool, dense gas at the centers of many 
galaxy clusters emits copiously at X-ray wavelengths. In the 
absence of a significant heat source, this gas would be ex- 
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pected to cool quickly and promote significant star forma- 
tion, at rates an order of magnitude larger than observed 
(see Peterson & Fabian 2006 for a review). The energy re- 
quired to offset this cooling is widely believed to come pre- 
dominantly from the central AGN, though the exact pro- 
cess or processes by which this occurs remains unclear (see 
e.g. McNamara & Nulscn 2007). AGN inflate bubbles of 
relativistic radio plEisma, which displace the X-ray emit- 
ting gas, forming clear cavities. This process is expected 
to provide a significant source of heating and turbulence 
(Churazov et al. 2001; see also e.g. Briiggen k, Kaiser 2002; 
Kaiser 2003; Briiggen 2003; De Young 2003; Ruszkowski et 
al. 2004a,b; Heinz & Churazov 2005). Sound waves and in- 
ternal waves generated by cavity inflation also provide a 
mechanism to heat the intra-cluster medium (ICM). These 
waves have been observed in both the Perseus and Cen- 
taurus clusters {e.g. Fabian et al. 2003; Fabian et al. 2006; 
Sanders & Fabian 2007; Sanders & Fabian 2008). AGN in- 
duced shock activity is a further promising avenue to heat 
the ICM. Clear indications of AGN-induced shocks are seen 
in several nearby galax;ies and clusters, the best examples 
being Hydra A, Virgo, and Perseus clusters {e.g. Nulsen 
et al. 2005a; Simionescu et al. 2009; Forman et al. 2005, 
2007; Fabian et al. 2003, 2006). AGN induced shocks in clus- 
ters are typically weak; Hercules A hosts the strongest ob- 
served AGN driven shock in a cluster with a Mach number 
M ~ 1.65 (Nulsen et al. 2005b). These outbursts can also be 
very energetic; MS 0735-1-7421 boasts the highest observed 
energy with ~ 6 x 10®^ ergs associated with its shock (Gitti 
et al. 2007). In galaxies and galaxy groups, these shocks can 
be much stronger with Mach numbers as high eis M ^ 8 as 
observed in Centaurus A (Croston et al. 2009; Kraft et al. 
2003). In M87, Forman et al. (2005) identified shock fronts 
associated with an AGN outburst 1 — 2 x lO'^ years ago. 
They argued that shocks may be the most significant mech- 
anism of heating the cool ICM near the core by the AGN. 
Forman et al. (2007) also measure temperature and den- 
sity jumps that are consistent with a weak shock with Mach 
number M ~ 1.2. However, important uncertainties remain. 
For example, significant temperature jumjjs are not seen 
across all putative shock fronts {e.g. Graham et al. 2008). 
Additionally, the observed pressure features are sometimes 
several kpc thick (Fabian et al. 2006), complicating our un- 
derstanding of how these processes work. 

In M87, two X-ray bright 'arms', discovered by Feigel- 
son et al. (1987) using the Einstein Observatory, are the 
most striking features driven by AGN activity. These arms 
extend to the southwest and east of the galactic center, are 
composed of cool gas, and are spatially coincident with ex- 
tended radio emission (Bohringer et al. 1995; Belsole et al. 
2001; Molendi 2002). Bohringer et al. (1995; see also Chura- 
zov et al. 2001) propose that these structures are the result 
of gas uplift in the wakes of bubbles of radio plasma buoy- 
antly rising in the hot ICM. However, the details of this 
process and the role of magnetic fields are largely not un- 
derstood. Subsequent Chandra and XMM-Newton observa- 
tions of M 87 have revealed many additional complexities in 
the system {e.g. Belsole et al. 2001; Matsushita et al. 2002; 
Molendi 2002; Bohringer et al. 2002; SakeUiou et al. 2002; 
Young et al. 2002; Forman et al. 2005, 2007; Simionescu et 
al. 2007, 2008; Werner et al. 2006). 

Here, we present the second in a series of papers study- 



Table 1. Summary of the Chandra observations of M87 used in 
this study. Columns hst the observation ID, detector, observation 
mode, exposure after cleaning and observation date. 



Ubs. lU 


Ubs. Date 


Detector 


Mode 


r!iXp. (ksj 


2707 


Jul. 6 2002 


ACIS-S 


FAINT 


82.9 


3717 


Jul. 5 2002 


ACIS-S 


FAINT 


11.1 


5826 


Mar. 3 2005 


ACIS-I 


VFAINT 


126.8 


5827 


May 5 2005 


ACIS-I 


VFAINT 


156.2 


5828 


Nov. 17 2005 


ACIS-I 


VFAINT 


33.0 


6186 


Jan. 31 2005 


ACIS-I 


VFAINT 


51.5 


7210 


Nov. 16 2005 


ACIS-I 


VFAINT 


30.7 


7211 


Nov. 16 2005 


ACIS-I 


VFAINT 


16.6 


7212 


Nov. 14 2005 


ACIS-I 


VFAINT 


65.2 



ing in unprecedented detail the thermodynamic properties 
of M 87. Using an ultra-deep (574 ks) observation, we uti- 
lize fully the excellent spatial resolution of the Chandra X- 
ray Observatory. This paper focuses on the overall thermo- 
dynamic structure and Fe distribution of the galaxy and 
surrounding cluster, and on the properties of the strongest 
shock features. Other papers will study the detailed physics 
of the X-ray bright arms and inner post-shock regions 
(Werner et al. 2010; hereafter Paper I) and the history of 
chemical enrichment (Million et al. 2010, in prep; Paper HI). 

The structure of this paper is as follows. Section 2 de- 
scribes the data reduction, basic imaging analysis, and spa- 
tially resolved spectroscopy. Section 3 presents detailed ther- 
modynamic maps for the cluster. Section 4 presents results 
on the radial properties of the ambient ICM to the north 
and south of the AGN; this excludes the X-ray bright arms, 
which introduce strong multi-temperature structure into the 
profiles. We place important new constraints on the isother- 
mality of the ambient cluster gas. Section 5 discusses the 
results, focusing on the shock features at radius r ~ 180 
arcsec, the isothermality of the ambient cluster gas, and a 
high significance metallicity bump at r 5 — 6 arcmin. Sec- 
tion 6 summarizes our conclusions. 

Throughout this paper, we assume that the cluster lies 
at a distance of 16.1 Mpc, for which the linear scale is 0.078 
kpc per arcsec. 



2 X-RAY OBSERVATIONS AND ANALYSIS 
2.1 Data Reduction 

Our work builds upon previous analyses of the same datasets 
studied here (see Young et al. 2002; Forman et al. 2005, 
2007). The Chandra observations were carried out using the 
Advanced CCD Imaging Spectrometer (ACIS) between July 
2002 and November 2005. The standard IeveI-1 event lists 
produced by the Chandra pipeline processing were repro- 
cessed using the CI AO (version 3.5.0) software package, in- 
cluding the appropriate gain maps and updated calibration 
products. Bad pixels were removed and standard grade se- 
lections applied. Where possible, the extra information avail- 
able in VFAINT mode was used to improve the rejection of 
cosmic ray events. The data were cleaned to remove periods 
of anomalously high background using the standard energy 
ranges recommended by the Chandra X-ray Center. The net 
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Figure 1. Exposure corrected, flat-fielded, surface brightness image of the central regions of the Virgo Cluster divided by the best-flt, 
radially symmetric doubIe-/3 model in the 0.6 — 2.0 keV energy band and smoothed with a 2 arcsec Gaussian. Unity represents locations 
where the data is equal to the model. This image uses only the ACIS-I data taken in 2005. The white circle represents the outer extent 
of the thermodynamic maps presented in Figs. [2l5l (radius 8.5 arcmin or ^ 40 kpc). 



exposure times after cleaning are summarized in Table [T] 
and total 574 ks. Separate photon-weighted response matri- 
ces and effective area files were constructed for each region 
analyzed. 



2.2 Imaging analysis 

Background subtracted images on a 0.984 x 0.984 arcsec^ 
pixel scale were created for each pointing in many narrow en- 
ergy bands (spanning 0.6 — 7.0 keV). These were flat-fielded 
with respect to the median energy for each image. Back- 
ground images were created from the blank-sky fields avail- 



able from the Chandra X-ray Center. The blank-sky fields 
were processed in an identical manner to the science ob- 
servations and were reprojected onto the same coordinates 
using the aspect solution files. The background images were 
normalized by the ratio of the observed and blank-sky count 
rates in the 9.5 — 12 keV band. The images were summed 
to form a series of combined background subtracted, expo- 
sure corrected images. Fig. [1] shows the X-ray image in the 
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Figure 2. Temperature map of the central regions of the Virgo Cluster (in units of keV) for the 8.5 arcmin (~40 kpc) radius region 
indicated by the white circle in Fig. [T] Regions have ~1,000 net counts, leading to la fractional uncertainties in the plotted quantities 
of ~10 per cent. Point sources, and readout errors have been excluded in this analysis and appear as black holes in the maps. 



0.6 — 2.0 keV band divided by the best-fit, double-/? modeQ 
and smoothed with a 2 arcsec filtered Gaussian. This im- 
age spans the entire region covered by ACIS-I observations 



taken in 2005. The white circle (r ^ 8.5 arcmin) denotes 
the outer extent of the data used for the spectral analysis 
reported here. 



^ In detail, this is the sum of two /3 models of the form 
\ 2" 



So 



1+^ 



, where So is the amplitude of the X-ray surface 



brightness, Tc is the core radius, and a is the power-law index. 



The cool X-ray arms are readily apparent in the image. 
An extensive imaging analysis of these data is presented by 
Forman et al. (2007). 
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Figure 3. Pressure map in units of keV cm X 2 Mpc j ' number density was calculated from the normalization of the spectra; 

rte = ^ , where K is the MEKAL normalization, A the projected area, and / the line-of-sight depth. C = 10^'* X 47rD£/1.2 and 
depends only on the distance to the source. We assume that the line-of-sight depth is i = 2 Mpc over the field of view. Other details as 
for Fig. [H 



2.3 Spatially- resolved spectroscopy 

2.3.1 Contour binning 

The individual regions for spectral fitting were determined 
using the contour binning method of Sanders (2006), which 
groups neighboring pixels of similar surface brightness until 



a desired signal-to-noise threshold is met. For data of the 
quality discussed here, regions outside the X-ray bright arms 
are small enough that the X-ray emission from each can 
be approximated usefully by a single temperature plasma 
model. Regions inside the X-ray bright arms show strong 
evidence for multi-phase gas (see Paper I). In this paper, 
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Figure 4. Entropy map of M 87 in units of kcV cm^ X ( 2Mpc ) ■ Other details as for Fig. |3] 



however, as a first approximation, we model all regions as 
a single temperature plasma, recognizing that this will bias 
the results for the X-ray bright arms. 

We constrain each region to have ~ 10'^ net counts, 
giving ~16,000 total regions for a total of ~16,000,000 net 
counts. This number of counts is sufflcient to determine the 
temperature of the cluster gas to < 10 per cent and the 
metallicity to ~20 per cent in each of the ~16,000 indepen- 
dent spatial bins. The density of the cluster gas is always 
determined to higher precision. We have imposed a lower 



limit on the region size of one arcsec'^, to mitigate the ef- 
fects of Chandra's point spread function. Regions vary in 
size from one to approximately 5,000 arcsec^ in the outer 
regions. Because the statistical uncertainty of the best-fit 
metallicity is ~ 20 per cent for regions with 1,000 counts, 
we have also performed a second analysis with regions con- 
taining at least 2,500 net counts each. This latter analysis 
produces statistical uncertainties of 3 — 5 per cent on the 
temperature and ~ 15 per cent on the metallicity. 
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Figure 5. Metallicity map in solar units. Abundances are with respect to the proto-solar values of Lodders et al. (2003). Regions contain 
at least 2, 500 net counts each. This results in statistical uncertainties of ~ 15 per cent per region. Regions inside the green contours 
contain significant 1 keV emission (see Fig. [2] Paper I) and are modelled improperly by a single temperature plasma. As a result, the 
metallicity in these regions is significantly underestimated (e.g. Buote et al. 2000; Simionescu et al. 2008; Werner et al. 2008). Other 
details as for Fig. [2] 



2.3.2 Annular binning 

To better study the properties of the ambient cluster gas, we 
have also performed a radial analysis in which we exclude 
the X-ray bright arms. This simplifies greatly the interpreta- 
tion of the results. For this analysis, we use annuli that vary 
in width from 10 to 30 arcsec, with ^250,000 net counts in 



each for a total of 36 annuli. Each annulus is fitted with a 
single temperature MEKAL or APEC spectral model (Sec- 
tion [2]3]3} with Galactic absorption. We have analyzed the 
northern (position angie^ -135 to 86 degrees) and southern 

^ All position angles are measured counter-clockwise from north. 
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(position angles 86 to 225 degrees) profiles both separately 
and together, always excluding the arms. 



2.3.3 Plasma models 

The spectra have been analyzed using XSPEC (version 12.5; 
Arnaud 1996). We use the MEKAL and APEC plasma emis- 
sion codes (Kaastra & Mewe 1993; Smith et al. 2001), and 
the photoelectric absorption models of Balucinska-Church & 



McCammon (1992). All spectral fits were carried out in the 
0.6 — 7.0 keV energy band. The extended C-statistic avail- 
able in XSPEC, which allows for background subtraction, 
was used for all fitting. 

The spectral model, applied to each spatial region, 
consists of an optically-thin, thermal plasma model, at 
the redshift of the cluster. We fix the Galactic absorp- 
tion to 1.93 X lO'^" atom cm~^ (determined from the 
Leiden/ Argentine/Bonn radio HI survey; Kalberla et al. 
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2005). The temperature, metallicity, and the normahzation 
are free parameters for every region (in detail, the normaliza- 
tions for each observation of each region are independently 
free). In determining the metallicity from the 2-dimensional 
spectral maps, all abundances are assumed to vary in a fixed 
ratio with respect to solar. In our higher signal-to-noise anal- 
ysis of annular spectra, (see Section [4|, we have allowed 
the abundances of O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni to 
be separate, free parameters. Results on the abundances of 
individual elements will be presented in a subsequent pa- 
per (Paper III). Abundances are given with respect to the 
'proto-solar' values of Lodders (2003) . 



2.3.4 Background modelling 

Background spectra for the appropriate detector regions 
were extracted from the blank-sky fields available from the 
Chandra X-ray Center. The blank-sky fields were processed 
in an identical manner to the science observation and were 
reprojected onto the same coordinates using the aspect so- 
lution files. Background regions were chosen to match the 
extraction regions for the science observations. The back- 
grounds were normalized by the ratio of the observed and 
blank-sky count rates in the 9.5 — 12 keV band (the statis- 
tical uncertainty in the observed 9.5 — 12 keV flux is less 
than 5 per cent in all cases). Due to the very small effective 
area of the Chandra telescope at hard X-ray energies, there 
is no significant cluster emission in the 9.5 — 12 keV band. 
Due to the high X-ray brightness of the target, uncertain- 
ties in the background modelling have little impact on the 
determination of thermodynamic quantities presented here. 



3 THERMODYNAMIC MAPS 

3.1 Large-scale thermodynamic properties 

Figs.[2HS]show the projected temperature, pressure, entropy, 
and metallicity maps determined from our analysis. Fig. [6] 
shows the same figures in mosaic form with the 90 cm (333 
MHz) radio contours from Owen et al. (2000) overlaid. Point 
sources and read-out errors have been masked and ignored. 
These results confirm spectrally, for the first time, many of 
the features identified by Forman et al. (2005, 2007) from 
an imaging analysis of the same data-set. 

The temperature map shown in Fig. [2] dramatically re- 
veals the complex temperature structure of the system. The 
X-ray bright arms stand out as relatively cool structures 
(Belsole et al. 2001; Young et al. 2002; Matsushita et al. 
2002; Molendi 2002; Forman et al. 2005; Simionescu et al. 
2007, 2008). They are spatially correlated with enhanced 90 
cm radio emission (Fig. [6^; contours are overlaid in white; 
from Owen et al. 2000). However, within ~3 arcmin (the 
radius of the shock front described in Section I5.2.1|l . the 
southwestern arm does not coincide precisely with the radio 
emission (see Fig. |6^; see also Forman et al. 2005, 2007) . 

Outside of the cool arms, the cluster appears remark- 
ably isothermal at a given radius (see also Section 5.1). 
The average temperature, excluding the arms, rises from 
~ 2 keV at a radius r ~ 0.9 arcmin (~ 4 kpc) to ~ 2.7 
keV by the edge of the field (r 8.5 arcmin). The X-ray 



bright, cold arms appear to be the most significant depar- 
tures from spherical symmetry. Surface brightness discon- 
tinuities at r ~ 3.75 (northwest of the AGN) and r ~ 8 
arcmin (south of the AGN) noted by Forman et al. (2005) 
do not have clear, corresponding features in these maps. 

The pressure map shown in Fig. [3] (Fig. [BJs with 90 cm 
radio contours from Owen et al. 2000) appears significantly 
more regular than the temperature map. The exception is 
a ring of high pressure at a radius of r 3 arcmin that 
is approximately ~ 0.7 arcmin thick. Forman et al. (2007) 
also see a feature at this radius in their analysis and argue 
that it is consistent with an AGN driven shock propagating 
through the ICM with a Mach number of M ~ 1.2. 

At the location of the cool arms, the observed pressure 
is somewhat lower. This is, however, likely an artefact due 
to the strong, projected multi-temperature structure in the 
direction of the X-ray bright arms in these regions (Paper 
I). 

The entropy map (Fig. |4]) reveals a similar morphologi- 
cal structure to the temperature map. The cool. X-ray bright 
arms have significantly lower entropy than the ambient ICM 
at the same radius. This suggests that the gas within the 
arms originates from the low entropy central regions of the 
cluster. The detailed thermodynamics of the arms are dis- 
cussed in Paper I. 

The metallicity of the ambient cluster gas (Fig. [5]) peaks 
at a value of ~ 1.5 solar near the core and drops to ~ 0.3 
near the edge of the field, r ~ 8.5 arcmin (40 kpc). The 
southern region is also more metal rich than the north and 
will be discussed further in Section 4.2. 

At first glance, our maps suggest a reduced metallicity 
within the X-ray bright arms with respect to the ambient 
gas (the green contours of 1 keV emission; see Paper I). 
However, this apparent decrement is likely a bias due to the 
complex multi-temperature structure within the arms {e.g. 
Buote 2000; Simionescu et al. 2008; Werner et al. 2008 and 
references within). The metallicity structure of the X-ray 
bright arms will be discussed further in Paper III. 

A high metallicity structure is observed at a radius of 
r ^ 4.5 arcmin to the west of the central AGN. This ridge 
connects with the high metallicity extension to the north- 
west of the core at a radius of r ~ 3.3 arcmin. However, this 
ridge is not continuous, with some regions of low abundance 
mixed within. As discussed in Section [5.31 this structure is 
likely due to the uplift of cool, metal-rich gas in the wake 
of buoyantly rising, radio-emitting plasma. High metallicity 
ridges are also observed in the Perseus Cluster, Ophiuchus 
Cluster, and Abell 2204 (see Sanders et al. 2005; Million et 
al. 2010; Sanders et al. 2009). 



3.2 Thermodynamic properties of the inner 
regions 

Temperature and pressure maps for the central 2x2 arcmin^ 
(9.4 X 9.4 kpc^) region of M87 are shown on an enlarged 
scale in Fig. [T] The 6 cm (5 GHz) radio contours of Hines 
et al. (1989) have been overlaid. Fig. [8] shows the 0.6 — 2.0 
keV residual X-ray image (extracted from Fig. [TJ for the 
same central 2x2 arcmin^ region. A number of prominent 
features in the images have been labelled with arrows. Each 
labelled cavity is detected at more than the 10a statistical 
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Figure 7. Thermodynamic maps (temperature on left; pressure on right) of the central 2x2 arcmin^ section of M 87. Radio contours 
from Hines et al. (1989) are overlaid in green. Several notable features are labelled with arrows. 

significance in surface brightness with respect to the smooth 
double-/? model. 

The jets, which have a hard non-thermal spectral en- 
ergy distribution (see e.g. Marshall et al. 2002; Wilson & 
Yang 2002; Perlman & Wilson 2005), would normally ap- 
pear as bright, high pressure, hot structures in the maps, 
but have been masked in these images. At a radius of 7 arc- 
sec (0.5 kpc) in the counter-jet direction, we observe a low 
X-ray surface brightness, hot cavity ('A') (see also Young 
et al. 2002). Interestingly, there is no obvious corresponding 
radio feature. A region of brighter, cooler gas encircles this 
structure and the AGN. 

Several other low surface brightness cavities are visible 
in Fig. 8, in the jet and counter-jet directions (see also For- 
man et al. 2007). These correspond to regions of apparent 
low pressure and high temperature. Two of these cavities (la- 
belled 'B' and 'C; both are at a radius of r ~ 40 arcsec from 
the nucleus in the jet and counter-jet directions) are also 
filled with radio emitting plasma. Rings of cooler gas sur- 
round these cavities. The eastern-most cavity ('D'), which 
also has relatively low pressure, lies at the edge of the radio 
cocoon that surrounds the AGN. Cavity E, labelled 'bud' by 
Forman et al. (2005), is located r ~ 40 arcsec (3 kpc) to the 
south of the AGN and is associated with high temperature, 
low pressure gas and is filled with radio plasma. Surrounding 
the cavity is a rim of high pressure material that is possibly 
shocked gas. 

The southwestern edge of the bright radio emission cor- 
responds to a high pressure, high temperature. X-ray bright 
'ridge'. We identify this ridge as a shock front. This feature 
(labelled 'Shock A') may extend to surround most of the 
bright, inner radio cocoon; it is seen clearly to the northwest, 
though is less obvious to the east. The connecting 'bridge' 
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between the central low temperature region and the south- 
west arm, which is also cool, appears 'interrupted' in pro- 
jection by this shocked gas. To the northwest, for position 
angles of -63—43 degrees, at a radius of r 40 arcsec from 
the AGN, a coherent arc of relatively cool gas is observed. 
This cool gas is a multi-phase structure and is discussed in 
detail in Paper I. The outer edge of this arc at r ~ 50 arcsec 
(labelled 'Cool Edge') is also probably physically linked to 
the central region of cooler gas, but is apparently broken up 
by the shock front ('Shock B' and 'Shock C). 



4 AMBIENT GAS PROPERTIES 

4.1 Azimuthally averaged projected profiles 

Our thermodynamic mapping reveals, in exquisite detail, 
the complex structure of the bright X-ray arms (see also 
Paper I). Just as importantly, Chandra's excellent spatial 
resolution also allows us to mask these structures and bet- 
ter determine the properties of the ambient gas surrounding 
M87. In this manner, we can avoid the complicating effects of 
temperature substructure, which will bias temperature and 
metallicity results in cases where the wrong thermal model 
is assumed {e.g. Buote 2000; Simionescu et al. 2008; Werner 
et al. 2008). 

For this analysis, we use 36 annuli that vary in width 
from 10 to 30 arcsec, each with ~ 250, 000 net counts. Each 
region is fitted in the first case with a single temperature 
MEKAL spectral model, with Galactic absorption. We also 
repeat the analysis using the APEC code. The high signal- 
to-noise of our annular spectra allows us to constrain the 
abundances of O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni simulta- 
neously, in addition to the temperature and normalization 
of the cluster gas (19 free parameters; 17 free parameters at 
large radius). Here we present results only for Fe. Results 
on elemental abundances other than Fe will be presented in 
a subsequent paper (Paper III). Error bars are drawn at the 
AC = 1 (68 per cent) confidence level. 

Fig. [9] shows the projected profiles of the temperature 
and Fe abundance. Results are shown for both the APEC 
(filled circles) and MEKAL (open squares) plasma codes. Our 
profiles extend to the edge of the mosaic field of view, r ^ 8.5 
arcmin (40 kpc) . The inner radius for these profiles is defined 
by the 'Cool Edge' identified at r ~ 50 arcsec (4 kpc) in Fig. 

The temperature profile (Fig. [9^) steadily rises from 
~ 2.1 keV to ~ 2.5 keV for 50 < r < 240 arcsec. Beyond 
r ~ 4 arcmin and out to r ~ 8.5 arcmin, the temperature 
rises more slowly but again steadily, reaching kT ~ 2.7 keV 
by r ~ 8.5 arcmin. The statistical errors on the measured 
temperatures range from 0.3 to 0.8 per cent. The results 
for the two plasma codes exhibit slight differences, with the 
MEKAL temperatures being ~ 2 per cent lower than those 
for APEC. 

Careful inspection of Fig. [9^ reveals a jump in temper- 
ature at a radius of r 130 arcsec. This corresponds to the 
inner edge of the shock front described in Section 13.11 

The Fe abundance profile (Fig. [SJj) peaks at Z-pc > 1.4 
solar in the central regions. The profile then exhibits a 
steady decline to Z-pc ~ 0.7 solar by r ~ 350 arcsec. A 
significant enhancement, or 'bump', in the Fe abundance is 



seen for radii 350 < r < 400 arcsec (27 — 31 kpc). This is 
approximately the radius at which the bright X-ray arms 
terminate (see also Section [5. 3p . 

We have investigated the impact of our choice of plasma 
code on the measured Fe abundances, finding a systematic 
offset of ~ 2 per cent in the results for the MEKAL and 
APEC codes at lower temperatures. This offset rises to ~ 10 
per cent for the hottest gas observed at the edge of the field. 

4.2 Southern and northern profiles 

We have analyzed separately the properties of the ICM to 
the north and south of the X-ray bright armslf] Here, we have 
used the APEC plasma code only. Fig. [TO] shows the pro- 
jected temperature and Fe abundance profiles to the north 
(position angles -135 to 86 degrees; filled circles) and south 
(position angles 86 to 225 degrees; open squares). Our re- 
sults agree well with those of Simionescu et al. (2007) using 
XMM-Newton data. 

The northern gas temperature profile is approximately 
constant at kT ~ 2.1 keV for 50 < r < 100 arcsec, and 
then rises to ~ 2.7 keV by 300 arcsec. Beyond this, the 
temperature rises only slightly out to r ~ 500 arcsec. In 
contrast, the southern profile shows a stronger temperature 
gradient. The central 100 arcsec dips below 2 keV whilst at 
large radii r ~ 450 arcsec the temperature rises to kT ~ 2.7 
keV. 

The measured Fe abundance is consistently lower in 
the northern sector, presumably due to metal transport by 
sloshing (see Simionescu et al. 2010). For both the north 
and south, however, the Fe abundance peaks in the cen- 
ter and drops steadily out to a radius r ~ 400 arcsec. To 
the north, we observe a bump of enhanced Fe abundance 
at 330 < r < 400 arcsec. A less significant bump is also 
observed to the south. This radius corresponds to the outer 
radius of the X-ray bright arms. A second large bump ap- 
pears to begin at r ~ 500 arcsec (39 kpc) to the north, 
which corresponds to the outer radius of the overall radio 
halo and suggests a similar origin. However, our profiles do 
not extend far enough to unambiguously confirm this fea- 
ture. The southern profile shows no convincing evidence of 
Fe abundance enhancements at r > 400 arcsec. 



5 DISCUSSION 

5.1 Isothermality of the ambient cluster gas 

The X-ray bright arms are well known to have multi- 
temperature structure {e.g. Belsole et al. 2001; Matsushita 
et al. 2002; Molendi 2002; Simionescu et al. 2008; Paper I). 
Previous spectral analyses have suggested that the gas out- 
side of the X-ray bright arms is approximately isothermal at 
a given radius (Matsushita et al. 2002). Our data, however, 
allow us to study the spatial temperature structure in more 
detail. Such measurements have important implications for 
the physics of the ambient ICM, including e.g. the efficiency 
of thermal conduction and the assumption of hydrostatic 
equilibrium. 

^ The cool X-ray bright arms remain excluded in this analysis. 
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Figure 9. Azimuthally averaged, projected profiles of temperature (left; in keV) and Fe abundance (right) with respect to Lodders 
(2003) proto-solar abundances. Results are shown for both the MEKAL (open squares) and APEC (filled circles) plasma emission codes 
(Kaastra & Mewe 1993, Smith et al. 2001, respectively). A radius of 100 arcsec corresponds to a physical scale of 7.8 kpc. The X-ray 
bright arms are excluded from this analysis (see text). The fits for these spectra include the abundances of O, Ne, Mg, Si, S, Ar, Ca, Fe, 
and Ni as separate, free parameters in addition to the temperature and normalization of the cluster gas. 




Radius (arcsec) Radius (arcsec) 

Figure 10. Projected profiles of the temperature (left; in keV) and Fe abundance (right; with respect to Lodders (2003) proto-solar 
abundances) for the northern and southern sectors, separately. Here we use the APEC plasma emission code only. Results for the northern 
sector (position angles -135 to 86 degrees) are plotted as filled circles. Results for the southern sector (position angles 86 to 225 degrees) 
are shown as open squares. 



To determine the presence of multi-temperature struc- 
ture at a given radius, we have made histograms of our differ- 
ent measurements of kT in 0.3 arcmin wide annular regions 
(Fig. I12|l . The analysis was performed in eight annuh in the 
radial range of 90 < r < 235 arcsec. For this analysis, we 
use position angles of -81-42 degrees. The width of the an- 
nuli was chosen to be large enough that ~ 75 different bins 
from the thermodynamic maps are present in each. For this 



analysis, we use the temperature map created with at least 
2, 500 net counts per region. The spectral models fitted to 
the individual spectra have the abundance ratios fixed to 
the best-fit values determined from the analysis described 
in Section|311 (see Section 5.1.1; see also Paper III). 

We have also tested each individual region for spectral 
evidence of multi-temperature structure. A 4 temperature 
model (see Paper I) shows an absence of any significant ther- 
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mal emission at temperatures between 0.5 < kT < 1.5 keV 
at larger than the few per cent level (see also Fig. Ilip . 

Having determined the distribution of temperature 
measurements, we compare this distribution with a Gaus- 
sian model. The width of the Gaussian is set to be equal 
to the average statistical error per bin (the typical standard 
deviations are approximately 3-5 per cent). Errors in each 
bin were determined from an MCMC analysis having chain 
lengths of at least 10* samples, after correcting for burn-in. 
A good match between the temperature measurements and 
the Gaussian model would be consistent with isothermality. 
Fig. [13] shows the ratio of the observed width of the emis- 
sion measure distribution (labelled a) divided by the expec- 
tation, given the observed mean statistical error (labelled 
5). 

Only one out of 8 regions show evidence for a deviation 
from isothermality. This region is located at the tempera- 
ture discontinuity behind the shock front. The temperature 
distribution in the other seven regions is consistent with be- 
ing isothermal. Projection effects are not expected to have 
a significant influence on the analysis. 

The observed isothermality at a given radius implies 
either that the ambient cluster gas beyond the arms is rela- 
tively undisturbed by AGN uplift, or that there is efficient 
conduction within a spherical shell. If conduction is respon- 
sible for the observed isothermality at a given radius, then 
the likely orientation of the magnetic fields is perpendicu- 
lar to the radial direction, as expected under action of the 
heat-flux buoyancy instability (Quataert 2008; see also e.g. 
Bogdanovic et al. 2009; Parrish et al. 2009). 

Visual inspection of Fig. 5 also suggests that the metals 
in the ambient ICM are significantly clumped. A detailed 
analysis of the metal distribution and a discussion of the 
implications for turbulence and conduction in the ICM will 
be presented in future work. 
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Figure 11. 6.5 X 6.5 arcmin^ section of tlic temperature map 
made with regions that contain ~ 2, 500 net counts. The 8 partial 
annuli used in the analysis are overplotted in green. The average 
error per bin ^ 5 per cent in temperature. 



Table 2. Summary of the shock models for each of the 4 differ- 
ent surface brightness profiles. Columns list the position angles 
of included data, shock radius, and Mach number of the observed 
shock within each position angle range. Conservative uncertain- 
ties on the Mach number, which bracket the observed scatter in 
the surface brightness profiles, are estimated to be M ± 0.04. 



5.1.1 Temperature bias 

Fig. [2] shows the temperature (left) and metallicity (right) 
profiles measured from both the peaks of the emission mea- 
sure distributions for the individual regions (open squares) 
and from the integrated annular spectra presented in Sec- 
tion [4] (filled circles). The temperature and metallicity re- 
sults from the two techniques are systematically different 
by 10 and 30 per cent respectively. The magnitude of this 
offset decreases as a function of radius. 

From spectral simulations, we have determined that this 
bias is a combined modelling and signal-to-noise issue. Low 
temperature spectra (below ~ 3 keV) contain a strong line 
contribution at energies where Chandra 's effective area is the 
largest. However, in low signal-to-noise spectra, one cannot 
determine the abundances of individual elements separately. 
These are instead implicitly tied to the Fe abundance. The 
flux from line emission is, therefore, modelled improperly 
which biases the temperature and metallicity to low values. 
Increasing the signal-to- noise of the spectra reduces these bi- 
ases, while simultaneously causing the residuals around the 
Fe L complex to significantly increase. The apparent posi- 
tion dependence is, therefore, explained by the temperature 
dependence of the intensity of these line components. The 
bias in low signal-to-noise spectra can be lessened by fitting 
a model with the abundance ratios of the individual de- 



position Angle (degrees) r (arcsec) M 



20-50 182 1.16 

190-220 168 1.29 

320-350 188 1.21 

350-20 180 1.34 

ments adjusted to their best-fit values, as determined from 
a fit to the higher signal-to-noise spectra. 

We emphasize that this bias does not affect the primary 
conclusions of this paper. Where appropriate, we have used 
the abundance ratios measured from the high signal-to-noise 
radial analysis (see Paper III) to reduce the impact of this 
effect. 

5.2 Analysis of the shock fronts 

5.2.1 The pressure ring at r ^ 180 arcsec (14 kpc) 

The large-scale pressure map (Fig. [3]) reveals the presence 
of a thick (~ 40 arcsec) high pressure ring at ~180 arcsec. 
This feature has previously been identified as a shock with 
Mach number M ~ 1.2 and age ~ 14 Myrs (Forman et al. 
2007). 

Fig. [TS] shows the projected surface brightness (upper 
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Figure 12. Histograms of temperature measurements (in iceV) for all 8 regions. Each region contains ~ 75 independent temperature 
measurements. Gaussians witli standard deviation equal to that of the average measurement error per bin are overplotted (dashed lines). 
Only one region shows significant deviation away from isothermality. 
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Figure 13. The ratio of the standard deviation of a Gaussian 
model fitted to the observed temperature distribution (cr) and the 
average error per bin in the thermodynamic maps (5) as discussed 
in Section 15.11 The ratio of errors is consistent with a constant 
value, which is overplotted on the data. 



panels; in cts s ^ arcsec ^) for the position angles 190- 
220 (south of the AGN; left) and position angles 350-20 
(north of the AGN; right). These surface brightness pro- 
files show signatures commensurate with a weak shock at 
a radius of r ~ 180 arcsec. The relevant features occur at 
a slightly smaller radius to the south with respect to the 
north (r ~ 170 as opposed to r ~ 180 arcsec). Overplot- 



ted on these profiles are simple shock models for a range of 
Mach numbers: for the northern sector, Mach numbers of 
1.30, 1.34 and 1.38; for the southern sector, Mach numbers 
1.24, 1.29 and 1.33. 

The shock models assume a spherically symmetric, 
hydrodynamic model of a point explosion in an initially 
isothermal, hydrostatic atmosphere. The initial gas density 
profile is assumed to be a power law, which is adjusted to 
match the observed surface brightness profile beyond the 
shock (Nulsen et al. 2005a; Simionescu et al. 2009). We have 
also explored a model that injects a constant amount of en- 
ergy per unit time into the system. This model provides no 
significant change to the calculated surface brightness pro- 
file, which is already fit well by the simple point explosion 
model. 

These models show good overall agreement with the 
data and the analysis performed by Forman et al. (2007). 
Table 2 shows the Mach number and shock radius from 4 
different 30 degree wedges (from -40-50 degrees and also 
from 190-220 degrees). Systematic uncertainties dominate 
the modelling of the surface brightness in these wedges. 
These include any projected non-sphericity of the shock due 
to the likely jet axis alignment near our line of sight, and the 
incorrect assumption of an initially isothermal, hydrostatic 
atmosphere. We, therefore, do not include standard, statis- 
tical error bars from the models. Instead, we estimate the 
uncertainty in the Mach number by finding a Mach number 
which brackets the observed scatter in the surface brightness 
profiles (see Fig. ll5|) . This represents a conservative estimate 
of the uncertainty in the Mach number of the shock. The 
overall shape of the surface brightness constrains the Mach 
number of the shock to within ~ 5 per cent. The radius and 
Mach number of the shock vary at the ten per cent level as a 
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Figure 14. Temperature (left) and metallicity (right) measured from the mean of the emission measure distribution when fit with a 
Gaussian (open squares) and from the integrated spectrum using 250,000 net counts and the same, simple spectral model for the same 
annular regions in Section |4] (filled circles). The spectral model is a single phase temperature model with elemental abundances tied 
to Fe. Both the temperature and metallicity are systematically offset by 10 and 30 per cent respectively when measured with the two 
techniques. This is due to the incorrect modelling of the line emission around the Fe L complex, where Chandra's effective area is largest. 



function of position angle. This is likely due to the detailed 
structure of the ambient gas and shock front in each sector. 

The temperature profiles for the position angles of 
— 135 — 86 and 86 — 225, drawn from Fig. 1101 are shown in 
the lower panels of Fig. [TS] We have approximated the am- 
bient temperature profile as a linear fit to the plotted data 
range, excluding the bins within the ~ 40 arcsec thick shock. 
Under this assumption, we observe that both the northern 
and southern projected temperature profiles contain at least 
four bins (between 130 ^ r < 170 arcsec) that appear hotter 
(~ 5 per cent) than expected. Fig. [TS] shows the expected 
projected temperature profile divided by the pre-shock tem- 
perature profile as a function of radius. The ^^5 per cent 
rise in the temperature and the ~40 arcsec thick shock in 
the observed, projected temperature profile are fully consis- 
tent with the shock model given an average Mach number 
M = 1.25 (see Fig. [T6)) . We, therefore, conclude that the 
high pressure ring is consistent with a shock with M = 1.25. 

Thick shocks have also been observed in the Perseus 
Cluster (see Fabian et al. 2006) , which does not show a clear 
temperature jump at the leading edge of the front. The small 
increase in temperature within the thick over-pressurized 
ring in M 87 suggests that shock heating is present in these 
systems but that the thickness and the small amount of 
heating associated with these weak shocks makes tempera- 
ture increases difficult to detect. 



5.2.2 The inner shock 

A second shock is seen closer to the core (~ 40 arcsec; Sec- 
tion 3.2; Fig. [7| The properties of this shock are difficult to 
constrain due to its proximity to the core. Many other sur- 
face brightness features {e.g. the cavities filled with radio 
plasma seen in Fig. [T]) are observed at similar radii, prevent- 



ing an accurate determination of the pre-shock gas proper- 
ties and the Mach number. However, the projected temper- 
ature jump seen to the south of the AGN at r ~ 40 arcsec (3 
kpc; Fig. O, suggests a Mach number of at least M > 1.2. 

The presence of the two shocks suggests that AGN out- 
bursts are fairly common. The estimated Mach numbers 
from these shocks and their relative distances can be used 
to determine the frequency of these outbursts. These shocks 
suggest that AGN outbursts occur approximately every ~ 10 
Myrs. Shock heating is likely only relevant in the central re- 
gions of the cluster because shocks significantly weaken as 
they expand. If repeated shocks occur every ~ 10 Myrs, as 
is suggested here, then AGN driven, weak shocks could pro- 
duce enough energy to offset the radiative cooling of the 
ICM (Nulsen et al. 2007). 



5.3 Uplift of cool, metal-rich gas 

Bohringer et al. (1995; see also Churazov et al. 2001) pro- 
pose that rising bubbles filled with radio emitting plasma 
may be responsible for dragging cool, metal-rich gas up out 
of the central regions of clusters. Buoyant bubbles will rise 
to a height that depends on their entropy, stopping when 
they reach their appropriate adiabat. This process may ex- 
plain the observed radio edges at r ~ 300 — 400 arcsec (23-31 
kpc) and also at r ~ 400 — 500 arcsec (31-39 kpc). When 
the bubbles reach their maximum height, they flatten. A 
consequence of this interpretation is that metal-rich gas ris- 
ing with the radio bubbles should be deposited at a radius 
similar to that at which the bubbles flatten. Our abundance 
profiles show evidence for enhancements in Fe abundance at 
320 < r < 400 arcsec (see Figs. [9b. [robl. This is a similar 
radius to the outer edges of the X-ray arms. Indeed, bumps 
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Figure 15. Projected surface brightness (top; in cts arcsec"'^) for position angles 190-220 (left) and position angles 350—20 (right). 
The projected temperature profiles (bottom; in keV) are drawn from Fig. 1101 (left is south; right is north). The surface brightness profiles 
show a clear discontinuity at ^ 170 — 180 arcsec, corresponding to the outer edge of the shock. For the southern sector (position angles 
190-220; left) shock models with Mach numbers 1.25, 1.29, and 1.33 at a radius of 168 arcsec are plotted; for the northern sector (position 
angles 350-20; right), models with Mach numbers 1.30, 1.34, and 1.38 at a radius of 180 arcsec are shown. The overplotted spread in 
Mach numbers represent a conservative estimate for the uncertainty in the measured Mach number. For the temperature profiles, we 
approximate the ambient temperature profile as linear fits to the temperature profiles in the plotted radial range, excluding temperature 
bins within the ~ 40 arcsec thick shock. The projected temperature in the region of the shock is ~5 per cent hotter than the best fit 
ambient model. 



are observed in the abundance profiles of all elements except 
Oxygen (Paper III). 

The amount of Fe that is uplifted can be determined 
from the excess metallicity measured at that radius and the 
gas mass. In detail, 

Mpe = Mg X AZpc X /fc, (1) 

where Mfc and AZp^ are the mass and excess abundance 



of Fe, Mg is the gas mass, and /fo = (mFo) / {rriiQ) x 
(iVpc/EiiVi)iodd is the mass fraction of Fe (Fodders 2003). We 
measure the gas mass and excess Fe abundance in the north- 
ern sector from a deprojection analysis (using data from po- 
sition angles of -108-55 degrees). This implies a total mass 
of uplifted Fe of ~ 1.0 x lO'' Mq. We compare this to the Fe 
mass currently present in the arms (Simionescu et al. 2008). 
For a gas mass of 5 x 10® Mq and an Fe abundance of 2.2, 
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those authors measure an Fe mass of 1.5 X 10® Mq in the 
X-ray bright, cool arms. Thus, a single generation of up- 
lifted metals from buoyant bubbles is, in principle, enough 
to explain the observed excess Fc mass at large radius. It 
is unknown how much of the uplifted metals will remain 
at large radius after the bubbles flatten. Lower significance 
bumps at higher radius that coincide with previous gener- 
ations of radio bubbles may argue that this process occurs 
over multiple generations. It is likely that the magnetic field 
configuration and its evolution will also have an impact {e.g. 
Bogdanovic et al. 2009). In the Perseus Cluster, metallicity 
enhancements also correspond well to the edges of several ra- 
dio lobes, which suggests a similar origin to that discussed 
here (see Sanders et al. 2005). This process may also partly 
explain the metallicity ridges observed around the central re- 
gions of some cool core galaxy clusters (Sanders et al. 2005; 
Sanders et al. 2009; Million et al. 2010). 

We note that bumps in the Fe abundance profiles are 
not clearly observed to the south. However, the mean Fe 
abundance is larger in the south than in the north and a 
cold front is located at a similar radius (Simionoscu et al. 
2007, 2010), making the identification of such a bump more 
difficult. Our results support the rising and pancaking bub- 
ble simulations presented in Churazov et al. (2001; see also 
e.g. Briiggen & Kaiser 2002; Kaiser 2003; Briiggen 2003; 
Ruzkowski et al. 2004a,b; Heinz & Churazov 2005). 



6 SUMMARY 

Utilizing a detailed, spatially-resolved spectral mapping and 
an ultra-deep (574 ks) Chandra observation of M 87 and the 
central regions of the Virgo Cluster, we present an unprece- 
dented, close-up view of AGN feedback in action. Our maps 
reveal X-ray bright arms that have been lifted up by buoyant 
radio bubbles as relatively cool, low entropy features that 
are rich in structure (see e.g. Belsole et al. 2001; Molendi 
2002). The detailed properties of these arms are discussed 
in Paper I. Outside the arms, the ambient cluster gas is 
strikingly isothermal at each radius. This suggests that ei- 
ther the gas remains relatively undisturbed by AGN uplift 
or that conduction is efficient along the azimuthal direction, 
as expected under action of the heat-flux buoyancy instabil- 
ity (HBI). Many cavities axe seen in the inner radio cocoon 
(~ 40 arcsec or ~ 3 kpc) including one located at ~ 7 arc- 
sec (~ 0.5 kpc) from the central AGN in the counter-jet 
direction. These cavities are associated with decreased sur- 
face brightness, high temperatures, low thermal pressure, 
and typically appear filled with radio plasma in projection. 

Our pressure map confirms spectrally previous results 
on the presence of a circular shock front at r ~ 180 arcsec 
(14 kpc). Based on the surface brightness jump, we mea- 
sure an average Mach number of M = 1.25 for the shock at 
r ~ 180 arcsec, in good agreement with the analysis of For- 
man et al. (2007). Under the assumption of a linearly rising 
temperature profile, we also measure a temperature jump 
within the ~ 40 arcsec thick shock, consistent with the ex- 
pectations for a shock with Mach number M — 1.25. The 
thickness and overall weakness of these shocks will make it 
difficult to detect temperature jumps unambiguously in hot- 
ter clusters like Perseus {e.g. Fabian et al. 2006; Graham et 
al. 2008). 
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Figure 16. Expected emission weighted projected temperature 
profile divided by tiie pre-shock temperature given tlie Mach num- 
bers determined from the shock models in 4 angular sectors. The 
average shock with Mach number M = 1.25 at a radius of r = 180 
arcsec is represented by the solid line. The observed ~5 per cent 
increase in temperature and the ~ 40 arcsec thickness of the over- 
pressurized region is well explained by the shock model. 



Another shock, associated with the central radio co- 
coon, is observed at ^ 40 arcsec (~ 3 kpc). Because of the 
complexity of this region, we can only estimate the Mach 
number of this shock, using the observed projected temper- 
ature jump across the front. This suggests a Mach number 
of M > 1.2. The relative distances and Mach numbers of 
these shocks suggest that AGN outbursts are common and 
occur every ~ 10 Myrs. 

Our results show that a significant mass of metals is 
lifted up to large radius in the wake of the radio bubbles. At 
least 1.0 X 1O*^M0 of Fe is uplifted to a radius of 350 - 400 
arcsec (27-31 kpc) to the north of the AGN. Given that a 
similar amount of Fe (1.5 x lO^M©) is present in the X-ray 
bright arms (Simionescu et al. 2008), a single generation of 
bubbles can, in principle, explain the observed Fe excess. It 
is still uncertain how much metal falls towards the central 
regions after this process occurs. Fe uplift of this magni- 
tude may explain the observed metallicity ridges observed 
around some cool core clusters. Studies of such abundance 
features can further enhance our understanding of transport 
processes. This will be the subject of future work. 
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